Potential energy surfaces for cluster emitting nuclei 



o 
o 

(N 

Or 
<D 

VO 
(N 



> 

o 

&\ 
O 

in 
o 



o 

=5 



X 



Dorin N. Pocnaru, 1,2, Radu A. Gherghescu, 1, 2 and Walter Greiner 2 

1 Horia Hulubei National Institute of Physics and Nuclear Engineering, 
RO-077125 Bucharest-Maqurele, Romania 



^Frankfurt Institute for Advanced Studies, J. W. Goethe Universitat, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 
(Dated: February 9, 2008) 

Potential energy surfaces are calculated by using the most advanced asymmetric two-center 
shell model allowing to obtain shell and pairing corrections which are added to the Yukawa-plus- 
exponential model deformation energy. Shell effects are of crucial importance for experimental 
observation of spontaneous disintegration by heavy ion emission. Results for 222 Ra, 232 U, 236 Pu 
and 242 Cm illustrate the main ideas and show for the first time for a cluster emitter a potential 
barrier obtained by using the macroscopic-microscopic method. 

PACS numbers: 24.75.-K, 25.85.Ca, 27.90.+b 



INTRODUCTION 

Recently we performed a systematic analysis Q of the 
experimental results concerning heavy particle radioac- 
tivities showing that our predictions within the ana- 
lytical superasymmetric (ASAF) model (see Q and the 
references therein) have been confirmed and that the 
strong shell effects of the daughter 208 Pb were not fully 
exploited. In this way we could make suggestions for the 
candidates to be used in the future experiments. 

In the present work we take advantage of using the 
most advanced two center shell model to study 
the potential energy surfaces (PES) of cluster emitting 
nuclei showing deep valleys due to the doubly magic 
heavy fragments 208 Pb and 132 Sn. The Strutinsky's Q 
macroscopic-microscopic method is used. A particularly 
deep valley is that of 208 Pb which proved to be of practi- 
cal importance not only for the production of superheavy 
nuclei but also for experimental search of cluster decay 
modes. Even for alpha decay it is possible to see such a 
valley if the emitter is 212 Po or 106 Te. In the later case 
the heavy fragment 102 Sn plays the important role. The 
potential barrier shape for a heavy ion decay mode may 
be obtained by cutting the PES at a given value of the 
mass and charge asymmetry. In this way one can com- 
pare the difference between the macroscopic barrier and 
the total one with shell and pairing corrections taken into 
account, providing a further justification of the ASAF 
barrier shape. 



shape, the potential energy is defined as 
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where E Y = a 2 A 2 ^ 3 {l - 3x 2 + (1 + l/x)[2 + 3x(l + 
x)] cxp(— 2/x)}, E® = a c Z 2 A~ 1 / 3 are energies corre- 
sponding to spherical shape and a 2 = a s (l — kI 2 ), 
I = (N — Z)/A, x = a/Ro, R = r A 1 / 3 . The parame- 
ters cicjhi, a c = 3e 2 /(5ro), and ro are taken from Moller 
et al. 

The relative Yukawa and Coulomb energies By = 
EyjE Yl B c = E c /E® are functions of the nuclear 
shape. The dependence on the neutron and proton 
numbers is contained in E Y , in the fissility parameter 
X = E^/(2E Y ) and By. For a binary fragmentation 
with charge densities p\ e and p2e, one has [8j a relative 
energy 
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with axially-symmetric shape-dependent terms expressed 
by triple integrals 
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MACROSCOPIC ENERGY 
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In a binary fission process A Z — > Al Z\ + A2 Z 2 the 
phenomenological energy Ey+EM is calculated within u _ A ± 

Yukawa-plus-exponential model (Y+EM) 0, Q by tak- 
ing into account the difference between charge and mass 



in which by = -d 4 (r /2a 2 )a 2 R A/E Y , d = {z"-z')/2R 
is the nuclear semilength in units of Rq and 



asymmetry ig. g x - x' dy 2 

By requesting zero deformation energy for a spherical F± — y + yy± cos tp - —— (6) 
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Q Y = {[VP{^fP + 2a/R d) + 2a 2 /(R d) 2 } ■ 

exg>{-RoVPd/a) - 2a 2 /(R a d) 2 }/P 2 (7) 

F 2 is obtained from F\ by replacing dy 2 /dx with dy\/dx'. 
In the above equations P = y 2 +yj — 2yyi cosip+(x— x') 2 , 
w = tp/2ir, and x c is the position of separation plane be- 
tween fragments with -1, +1 intercepts on the symmetry 
axis (surface equation y = y(x) or y\ = y(x')). The 
integrals are computed numerically by Gauss-Legcndrc 
quadratures. 

In a similar way the Coulomb relative energy is given 

by 

B c = § u =(^) 2 B cl + ^B cl2+ ( P -^) 2 B c2 (8) 
and for axially symmetric shapes 

/x c px c 
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where b c = 5d 5 /8ir. In the integrand 

F(x,x')= {y yi [(K-2D)/3}- 
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K and K' are the complete elliptic integrals of the first 
and second kind, respectively: 



K(k) 



I (1 - k 2 smH)-^ 2 dt 
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K'(k) = (1 - k 2 sin 2 t) 1/2 dt 
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and a 2 p = (y + Vl ) 2 + (x - x') 2 , k 2 = 4yyi/a% D = 
(K — K')/k 2 . The elliptic integrals may be calculated by 
using Chebyshev polynomial approximation. For x = x' 
the function F is not determined. In this case, after 
removing the indetermination, we get F(x,x') = 4y 3 /3. 

Starting from the touching point configuration, R > 
i?t, for spherical shapes of the fragments, one can use an- 
alytical relationships. The Coulomb interaction energy of 



a system of two spherical nuclei, separated by a distance 
R between centers, is E c i2 — e 2 Z\Zi/ R, where e is the 
electron charge. 

Within a liquid drop model (LDM) there is no con- 
tribution of the surface energy to the interaction of the 
separated fragments; the barrier has a maximum at the 
touching point configuration. The proximity forces act- 
ing at small separation distances (within the range of 
strong interactions) give rise in the Y+EM to an inter- 
action term expressed as folllows 
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In many cases the interaction energy is maximum at a 
certain distance R m > Rt = R1 + R2, which can be found 
by solving numerically the following nonlinear equation 



e x + pi + x(pi + xp) = ; x = R/c 



in which 
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and the interval x t = Rt/a, x t + 5 may be given as in- 
put data of a program using Miillcr's iteration scheme of 
succesive bisections and inverse parabolic interpolation. 



SHELL AND PAIRING CORRECTIONS 

In the following we would like to outline the calcula- 
tions of the shell 0] and pairing 0| corrections SE = 
5U + 8P leading to the total deformation energy 

E de} (R, V ) = Ey+em (R, V) + SE{R 7 V ) (21) 

More details are given in the book By choosing two 
intersected spheres for nuclear shape parametrization one 
can take the separation distance between fragment cen- 
ters, R, as a deformation parameter. Initially, for a 
parent nucleus Ri = Rq — R 2 . At the touching point 
Rt = Ri + i?2. The mass asymmetry 77 = (A\ — A2)/A. 

The two-center shell model Q gives at every pair of 
coordinates (R, rj) the sequence of doubly degenerate 



3 



discrete energy levels e.; = Ei/hujQ in units of frio® = 
arranged in order of increasing energy. The 
smoothed-level distribution density is obtained by aver- 
aging the actual distribution over a finite energy interval 
r = jhujg, with 7—1, 

m = {El=i [2-1875 + W ( tt (1.75 - yi/6) 

-4.375)]e-»*} (1.772453857)" 1 (22) 

where y = x 2 = [(e — Ei)/j] 2 . The summation is per- 
formed up to the level n m fulfilling the condition \x%\ > 3. 
The Fermi energy, A, of this distribution is given by 



N p = 2 / g(e)de 



(23) 



with N p = Z for proton levels and N p = A — Z for neu- 
tron levels, leading to a non-linear equation in A, solved 
numerically. The total energy of the uniform level distri- 
bution 
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In units of Twjq the shell corrections are calculated for 



each pair (R, rf): 



Su(n, R, r,) = E 2ei(i2, i]) - u(n, R, rj) (25) 
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FIG. 1: Shell and pairing corrections for a symmetric (jj = 0) 
fission of 236 Pu. 

Similarly, for pairing corrections we take the doubly 
degenerate levels {e^} in units of TwJq. Z/2 levels are 
occupied with n levels below and n' above Fermi energy 
contributing to pairing, n = n' = Qg s /2. The cutoff 



energy, CI ~ 1 > A = 12/y/AHu>§. The gap A and Fermi 
energy A are solutions of the BCS system of two eqs: 

= Y , (k ~ - (26) 

V V(^-A) 2 + A2 
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where h = Z/2 - n + 1, k f = Z/2 + n', and 

As a consequence of the pairing correlation, the lev- 
els below the Fermi energy are only partially filled, while 
those above the Fermi energy are partially empty. Occu- 
pation probability by a quasiparticlc (uk) or hole (vk) is 
given by 

,,2 



u k - [1 - (c fc - \)/E k ] /2; ui = 
The quasiparticle energy is expressed as 
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(29) 



E u = v /(e,-A)2 + A2. (30) 

The pairing correction dp = p — p, represents the dif- 
ference between the pairing correlation energies for the 
discrete level distribution 

kf Z/2 A2 

(31) 

k=ki 

and for the continuous level distribution 



p = E 2v l ek ~ 2 E ek - ~a 
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p = -(gA 2 )/2 = -(g s A*)/4 



(32) 



Compared to shell correction, the pairing correction is 
out of phase and smaller (see FigQ leading for r) = con- 
stant to a smoother total curve Se(R) = 5u(R) + 5p(R) 
where Sp = Sp p + 8p n . 



RESULTS 

In the following wc shall present results for 222 Ra, 232 U, 
236 Pu, and 242 Cm, which are emitters of 14 C, 24 Ne, 28 Mg, 
and 34 Si, respectively. All have been experimentally ob- 
served. In order to obtain a relatively smooth PES we 
made the approximation rj = if a = {A\ — Ai)/A ~ ifz = 
(Z\ — Z<i)jZ. In this way the fragment nucleon numbers 
N\,Z\ and N 2l Z 2 (plotted in figures El EH HSI EJ are 
linear functions of r\. We prefer to use the dimensionless 
separation distance £ = (R — R4)/ (R t — R4) instead of R. 
In this way one can clearly see the initial parent nucleus 
at £ = and the touching point configuration at £ = 1. 
We adopt the usual convention of having zero deforma- 
tion energy and shell plus pairing corrections for the ini- 
tial spherical shape, leading to Edef = Ey+em = SE = 
at R — Ri for all values of g and at g = ±1 for all values 
of R. 
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the sixtieth by Walter Greiner to motivate the search for 
superheavy elements, and the development of Heavy Ion 
Physics worldwide and in Germany, where GSI was built. 
These valleys may be also seen on the total PES at the 
top of Fig. |2 Here the deepest valley remains that at a 
small value of rj not far from the minimum of the macro- 
scopic Y+EM energy at r) = 0, which is responsible for 
the cold fission. At a large value of 77, the 208 Pb + 14 C 
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FIG. 3: Contour plot of shell and pairing corrections for Ra 
vs (R - Ri)/(R t - Ri) and r). 



FIG. 2: PESof 222 Ravs(7?-7? i )/(^t-^) > and 77 = (Ai - 
A 2 )/(Ai+A2). Y+EM (bottom), Shell + Pairing corrections 
(center), and total deformation energy (top). The energies 
are expressed in MeV. 



The PES versus the normalized separation distance 
(R — Ri)/ (Rt — Ri) and the mass asymmetry 77 = {A\ — 
+ A 2 ) are plotted in Fig. The macroscopic 
Y+EM deformation energy is shown at the bottom, fol- 
lowed by the microscopic shell plus pairing corrections 
(center), and their sum (the total deformation energy) 
at the top. 

Three valleys around 77 ~ 0.8; 0.3 and 0.1 can be seen 
in the center of Fig. [2and on the corresponding contour 
plot (Fig. I2J- We only count the number of valleys for 
?7 > because the mirror rj < gives the same num- 
ber for complimentary heavy fragments becoming light 
ones and vice- versa. In figure El we shall see that they 
are produced due to the magicity of the nucleon num- 
ber of the fragments. Such cold valleys were used in 
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FIG. 4: A cut through the PES of 222 Ra at symmetry 77 = 
(top) and for 14 C radioactivity with 208 Pb daughter (bottom). 

valley, favouring the 14 C radioactivity of 222 Ra, is laying 
on the Businaro-Gallone mountain, hence it is shallower 
despite the fact that it is very pronounced in the shell 
correction surface. On the contour plot (figure |3J one 
can see how it evolves from lower mass asymmetry at a 
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FIG. 5: Deformation energies at the touching point configu- 
rations (R = Rt) of 222 Ra vs the asymmetry r\: Ey+em and 
Edef (top); shell + pairing corrections and only shell correc- 
tions (bottom). 
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FIG. 6: Variation of the neutron and proton numbers of 
the two fragments with the mass asymmetry at the touch- 
ing point, R= Rt for 222 Ra. 



small value of separation distance, £, to the larger rj at 
the touching point, £ = 1. 

Two plots obtained by cutting the total PES and that 
of Y+EM at a given value of the asymmetry parameter 
are shown in Fig.0] In the upper part one can see a two 
humped barrier at ?y = 0. The typical example for 14 C 
emission from 222 Ra, shown at the bottom, provide justi- 
fication for one of the basic assumption of the analytical 
superasymmetric fission model, which was very successful 
in predicting the half-lives of cluster decay modes. It is 
remarkable to see for the first time for a cluster emitter 
a potential barrier obtained by using the macroscopic- 
microscopic method. Having a smaller height and width 
compared to the (dotted line) macroscopic Y+EM bar- 
rier, it is very similar to the barrier used in ASAF which 
was lower and narrower than the Myers-Swiatecki's 
LDM barrier. 

One can see how deep are the two main valleys on the 
PES by plotting in Fig. 0a cut of the PES at the touching 
point configuration, R = R t . In the upper part of this 
figure the macroscopic energy, Ey+ e (the smooth dashed 
line), and the total deformation energy, E (heavy line) 



can be seen. In the lower part of Fig.[S]there are two other 
curves representing the shell plus pairing correction (dot- 
ted curve) and only the shell corrections (dashed line). 
As we already mentioned when the Fig. ^ was discussed, 
there are two effects of pairing corrections leading to a 
smoother variation and to a shallower valley. 
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FIG. 7: PES of 232 U vs (R - Ri)/{R t - R*) and jj = (Ax - 
A 2 )/(Ax+A 2 ). Y+EM (bottom), Shell + Pairing corrections 
(center), and total deformation energy (top). The energies 
are expressed in MeV. 

We can better understand how the deep valleys ap- 
pear in figure [3] if we plot in figure [B] the variation of the 
neutron and proton numbers of the two fragments with 
the mass asymmetry at the touching point, R = R t , for 
222 Ra. Every time a nucleon number reaches a magic 
value, the corresponding shell correction has a local min- 
imum. The very deep valley at rj > 0.8 (for different mass 
and charge asymmetry at t]a = 0.874 and r\z = 0.864) is 
produced by reaching almost simultaneously three magic 
numbers N\ = 126, N 2 = 8, and Z\ = 82 (the decay 
222 Ra — > 208 Pb + 14 C). The next one, at an intermedi- 
ate value of rj (t]a = 0.369, r\z = 0.363), is mainly due 
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to the Ni (Z 2 = 28) light fragment ( 222 Ra -> 152 Nd + 
2§Ni42). Finally the valley at a small asymmetry param- 
eter (t]a = 0.189 and r] Z = 0.136) corresponds to the cold 
fission process ( 222 Ra — > 132 Sn + 90 Sr). 
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top. There are again three valleys on the shell plus pair- 
ing corrections and on the contour plot (Fig. |5J) around 
77 ~ 0.8; 0.3, and 0.15. Unlike on the figure El for 222 Ra 
in which the valley due to Sn at a low mass asymmetry 
is not so deep, now this cold fission valley is well shaped 
and the trend continues for heavier nuclei such as 236 Pu 
and 242 Cm. An intermediate valley is produced by the 
magic neutron number N2 = 50 of the light fragment. 
They are also present on the total PES at the top of 
Fig. Here the deepest valley remains that of Sn, which 
is responsible for the cold fission; the 208 Pb + 24 Ne val- 
ley, explaining the 24 Nc radioactivity of 232 U, is laying 
on the Businaro-Gallonc mountain, hence it is shallower. 



FIG. 8: Contour plot of shell and pairing corrections for 232 U 
vs (R - Ri)/(Rt - R4) and rj. 




0.0 0.5 1.0 1.5 



2.0 



> 20 

CD 

|10 




"D 




LU 




30 




20 




10 










0.0 0.5 1.0 1.5 2.0 

(R - Ri)/(R t - Ri) 

FIG. 9: A cut through the PES of 232 U at symmetry 77 = 
(top) and for 24 Ne radioactivity with 208 Pb daughter (bot- 
tom). 



FIG. 10: Deformation energies at the touching point config- 
urations (R = Rt) of 232 U vs the asymmetry 77: Ey+em and 
Edef (top); shell + pairing corrections and only shell correc- 
tions (bottom). 
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FIG. 11: Variation of the neutron and proton numbers of 
the two fragments with the mass asymmetry at the touching 
point, R = R t for 232 U. 



232 jj 

Three PES versus the normalized separation distance 
and the mass asymmetry are plotted in Fig. the 
macroscopic Y+EM deformation energy at the bottom; 
the microscopic shell plus pairing corrections at the cen- 
ter, and their sum (the total deformation energy) at the 



In the upper part of Fig. one can see a two humped 
barrier at 77 = 0. The example for 24 Ne emission from 
232 U, shown at the bottom, provide justification for one 
of the basic assumption of the analytical supcrasymmct- 
ric fission model already mentioned in the previous sub- 
section. 

A cut of the PES at the touching point configura- 
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tion, R — R t is shown in Fig. ^] with the Ey+ em (the 
smooth dashed line) and the total deformation energy 
Edef (heavy line) in the upper part. The three valleys 
mentioned above are present both on the total deforma- 
tion energy (top) and in the lower part where the shell 
plus pairing correction (dotted curve) and only the shell 
corrections (dashed line) are plotted. 




mass and charge asymmetry (j]a = 0.138, r\z = 0.087). 
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FIG. 12: PES of 236 Puvs (R - Ri)/(R t - Ri) and r) = (Ai - 
A 2 )/(A 1 +A 2 ). Y+EM (bottom), Shell + Pairing corrections 
(center), and total deformation energy (top). The energies 
are expressed in MeV. 

The figure 1111 showing the variation of the neutron 
and proton numbers of the two fragments with the mass 
asymmetry at the touching point, R = R t , for 232 U al- 
lows us to understand how the deep valleys from figure lTUI 
are produced. The very deep valley around r\ = 0.8 
(i]A = 0.793, r]z = 0.783) is produced by the negative 
shell corrections due to two magic numbers Ni = 126 
and Z 1 = 82 ( 24 Ne radioactivity 232 U -> 208 Pb + 24 Ne). 
Another doubly magic heavy fragment, 132 Sn, produces 
the cold fission valley ( 232 U ->■ 132 Sn + 100 Mo) at low 



FIG. 14: A cut through the PES of 236 Pu at symmetry r\ = 
(top) and for 28 Mg radioactivity with 208 Pb daughter (bot- 
tom). 

The intermediate one (77^4 = 0.293, t]z = 0.304), quite 
shallow, is due to the neutron magic number of the light 
fragment N 2 = 50 ( 232 U -> 150 Nd + l|Ge 50 ). 

236p u 

PES versus the normalized separation distance and the 
mass asymmetry are plotted in Fig. El where Y+EM 
deformation energy is shown at the bottom, the shell 
plus pairing corrections at the center, and their sum (the 
total deformation energy) at the top. Very deep valleys 
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due to the doubly magic fragments 208 Pb and 132 Sn can 
be seen in the center of Fig. 1 121 and on the corresponding 
contour plot (Fig. 1 1 3|> around rj ~ 0.75 and 77 < 0.1, 
respectively. At an intermediate value of the asymmetry 
parameter, rj ~ 0.28 there is another valley not so deep. 
These valleys may be also seen on the total PES at the 
top of Fig. ^] Here the deepest valley remains that of 
Sn, which is responsible for the cold fission; the 208 Pb + 
28 Mg valley, explaining the 28 Mg radioactivity of 236 Pu, 
is laying on the Businaro-Gallone mountain, hence it is 
shallower. 

Two plots obtained by cutting the PES at a given value 
of the asymmetry are shown in Fig. 1141 a two humped 
barrier at r\ = in the upper part and the barrier for 2 Ne 
spontaneous emission from 236 Pu, shown at the bottom. 
The last one provides a qualitative microscopic justifica- 
tion for the chosen barrier shape within the phenomeno- 
logical ASAF. 




FIG. 15: Deformation energies at the touching point configu- 
rations (R = Rt) of 236 Pu vs the asymmetry r\: Ey+em and 
Edef (top); shell + pairing corrections and only shell correc- 
tions (bottom). 
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FIG. 16: Variation of the neutron and proton numbers of 
the two fragments with the mass asymmetry at the touching 
point, 7? = Rt for 236 Pu. 

One can see how deep are the three valleys by plot- 
ting in Fig. ED a cut of the PES at the touching point 
configuration, R = R t . In the upper part of this fig- 



ure the macroscopic energy, Ey+E (the smooth dashed 
line), and the total deformation energy, Edef (heavy line) 
can be seen. They exhibit few valleys of which that of 
Sn at lower values of 77 and that of Pb at higher 77 are 
the deepest ones. In the lower part of this figure there 
are two other curves representing the shell plus pairing 
correction (dotted curve), SE s h+ p , and only the shell cor- 
rections (dashed line), 5E s heii- The valley at the inter- 
mediate value of the asymmetry parameter, 77 ~ 0.28, is 
quite shallow. 
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FIG. 17: PES of 242 Cm vs (R - Ri)/(Rt - ft) and n = {Ai- 
A 2 )/(A 1 +A 2 ). Y+EM (bottom), Shell + Pairing corrections 
(center), and total deformation energy (top). The energies 
are expressed in MeV. 

The variation of the neutron and proton numbers of the 
two fragments with the mass asymmetry at the touching 
point, R = R t , for 236 Pu (figure lPo|) shows the occurence 
of magic numbers which generate the minima of the shell 
and pairing corrections producing the valleys on PES. 
The very deep valley around 77 = 0.75 (tja = 0.763, rjz = 
0.745) is produced by reaching simultaneously two magic 
numbers iVj = 126 and Z x = 82, 236 Pu -► 208 Pb + 
28 Mg. The cold fission valley ( 236 Pu -> 132 Sn + 104 Ru) 
at a small asymmetry (tja = 0.119, i]z = 0.064) is the 
result of two magic numbers Z\ = 50 and N\ = 82. The 
intermediate valley (i]a = 0.288, r\z — 0.277) is produced 
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FIG. 18: Contour plot ol shell and pairing corrections for 
242 Cm vs (R - Ri)/{R t - Ri) and rj. 

by the magicity of the light fragment N% = 50, 236 Pu 
-> 152 Nd + |*Se 5 o. 



There are four valleys on the PES of figureElfor 242 Cm 
at -q ~ 0.7; 0.4; 0.3, and 0.1. They are also present on 
the corresponding contour plot (Fig. I18|) . At the top of 
Fig. El the deepest valley remains that of cold fission 
due to a doubly magic heavy fragment 132 Sn, and the 



0.0 0.5 1.0 1.5 




0.0 0.5 1.0 1.5 

(R - Ri)/(R t - R^ 

FIG. 19: A cut through the PES of 242 Cm at symmetry r\ = 
(top) and for 34 Si radioactivity with 208 Pb daughter (bottom). 



silicon-and-lead valley the 208 Pb + 34 Si, explaining the 
34 Si radioactivity of 242 Cm. 




-0.8 -0.4 0.0 0.4 0.8 

V 



FIG. 20: Deformation energies at the touching point configu- 
rations (R = Rt) of 242 Cm vs the asymmetry rj: Ey+em and 
Edef (top); shell + pairing corrections and only shell correc- 
tions (bottom). 




5 40 
Z 20 



0.0 0.2 0.4 0.6 0.8 1.0 

n= (A, - A 2 )/A 

FIG. 21: Variation of the neutron and proton numbers of 
the two fragments with the mass asymmetry at the touching 
point, R = Rt for 242 Cm. Every time a nucleon number 
reaches a magic value, the corresponding shell correction has 
a local minimum. The very deep valley around 77 = 0.719 
is produced by reaching simultaneously three magic numbers 
Ni = 126, N 2 = 20, and Z x = 82. 

Two plots obtained by cutting the PES at a given value 
of the asymmetry are shown in Fig. 1191 In the upper 
part one can see the two humped barrier and at the bot- 
tom the barrier for 34 Si emission from 242 Cm. Having a 
smaller height and width compared to the (dotted line) 
macroscopic Y+EM barrier, it is very similar to the bar- 
rier used in ASAF which was lower and narrower than 
the the Myers-Swiatecki's LDM barrier. 

One can see how deep are the four valleys on the PES 
by plotting in Fig. |20l a cut of the PES at the touch- 
ing point configuration, R = R t . In the upper part of 
this figure the macroscopic energy, Ey+E, (the smooth 
dashed line) and the total deformation energy, E, (heavy 
line) can be seen. In the lower part there are two other 
curves representing the shell plus pairing correction (dot- 
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ted curve), SE s h+ p , and only the shell corrections (dashed 
line), SEaheii- 

The figure 1211 shows the variation of the neutron and 
proton numbers of the two fragments with the mass 
asymmetry at the touching point, R = R t , for 242 Cm. 
Like in the preceding subsections, the deepest valley on 
the shell and pairing corrections is produced by the very 
strong shell effect of the heavy fragment 208 Pb ( 242 Cm 
-> 208 Pb + ?|Si 20 ) at rj A = 0.719 and j] Z = 0.708; 
here there is also a contribution coming from the neu- 
tron number of the light fragment. N 2 = 20. The deepest 
valley on the total deformation energy is the cold fission 
one at r) A = 0.091, i] Z = 0.042, where 242 Cm -> 132 Sn + 
110 Pd. The two intermediate valleys are produced mainly 
by N 2 = 50 and Z 2 = 28. One has at r) A = 0.306, 
r)z = 0.292 the decay 242 Cm -> 158 Sm + ||Se 50 and at 
r] A = 0.421, m = 0.4 1 7 242 Cm -> 172 Er + ^Ni 42 . 

In conclusion the strong shell effect associated to the 
doubly magic character of the daughter 208 Pb, observed 
in the systematic analysis of experimental results 0, 
comes from a valley present on the potential energy sur- 
faces of cluster emitters at a relatively high value of the 
asymmetry parameter rj ~ 0.7 — 0.8. Despite its high 
depth on the microscopic corrections PES, on the total 
deformation energy it appears shallower since it is added 
to the Businaro-Gallone mountain of the Y+EM macro- 
scopic energy. The potential barrier shape of heavy ion 
radioactivity obtained from the first time by using the 
macroscopic-microscopic method provides further sup- 
port for the particular choice of the barrier within the 
analytical superasymmetric model. The depth of the cold 
fission valley produced by the doubly magic heavy frag- 
ment 132 Sn, which is small for 222 Ra, increases with the 
mass number of the parent nucleus being comparable to 



the lead valley on the microscopic corrections PES of 
242 Cm. Some other intermediate shallower valleys are 
produced by the magicity of the light fragment: Z 2 = 28 
for 222 Ra parent; N 2 = 50 for 232 U parent; N 2 = 50 for 
236 Pu; N 2 = 50 and Z 2 = 28 for 242 Cm. 
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